Temperature-dependent Raman scattering of natural and isotopically substituted PbS 
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Lead sulfide is an important semiconductor that has found technological applications for over a 
century. Raman spectroscopy, a standard tool for the investigation and characterization of semicon- 
ductors, has limited application to this material because of the forbidden nature of its first order 
scattering and its opacity to visible lasers. Nevertheless, useful vibrational spectra from two-phonon 
processes are obtained with red lasers, probably because of a resonance in the concomitant elec- 
tronic transitions. Here we report temperature dependent spectra, covering the 10-300 K range, for 
two samples with different sulfur isotopic compositions. The results are analyzed by comparison 
with ab initio calculations of the lattice dynamics of PbS and the corresponding densities of one 
and two-phonon states. Emphasis is placed on the analysis of the two phonon band centered at 
~ 430 cm -1 . 



PACS numbers: 

I. INTRODUCTION 

Lead sulfide (PbS), the mineral galena, is an impor- 
tant material in lead smelting [1. It can be found in 
nature and also grown as an n- or p-type semiconduc- 
tor with carrier concentrations as low as 10 17 cm -3 . Its 
small electronic gap (0.4 eV), with an anomalous temper- 
ature dependence [5] , implies small effective masses which 
make PbS useful as a material for nanostructures[3]. On 
the negative side, it is a visually offensive black product 
which appears during the degradation of lead white pig- 
ment in artwork It has been used as a nanoparticle 
source in very ancient metallic lustre pottery, possibly 
the oldest known form of controlled nanotechnology |5J. 
It has also found an age old application in the adulter- 
ation of an eye-lining cosmetic called "kohl" (Sb2S3)[B]. 

Raman spectroscopy provides a versatile, non- 
destructive technique for the investigation and character- 
ization of semiconducting materials [7J. The NaCl struc- 
ture of PbS, however, makes the material first order Ra- 
man forbidden. Nevertheless, it has been recently dis- 
covered that characteristic Raman spectra are obtained 
when lasers in the red are used as Raman spectroscopic 
sources. These spectra are dominated by a band (ac- 
tually a doublet) centered at around ~ 430 cm -1 and a 
somewhat narrower peak at ~ 210 cm . A few other 



characteristic weak structures are also observed. These 
spectra have been recently used for the characterization 
of lustre pottery and possibly adulterated "kohl" E] . 

In a recent publication [8 , our group investigated the 
dependence of the two aforementioned Raman bands on 
the isotopic mass of sulfur and on temperature, the lat- 
ter being limited to the 100-300 K region. The lack of 
data at lower temperatures hampers the investigation 
of the mechanisms responsible for the observed anhar- 
monic temperature shifts. Therefore, in the present work 
we undertook the investigation of these Raman spec- 
tra down to a temperature of ~ 10 K. As in the case 
of the anomalous temperature dependence of the elec- 
tronic gap (but most likely unrelated to it) the frequency 
of the ~ 430 cm" 1 band anomalously increases with in- 
creasing temperatures, a fact that suggests interaction of 
each of the two phonons responsible for this ~ 430 cm -1 
band with a set of two phonons at lower frequencies 
(~ 65 cm -1 ). The dependence of the 430cm" 1 band on 
the isotopic mass of sulfur, however, is normal and sug- 
gests that the corresponding eigenvectors are sulfur dom- 
inated (for a general discussion of the dependence of an- 
harmonic self-energies on isotopic masses see Ref. 9). In 
order to complete the assignment of the observed Raman 
structures, we have calculated the phonon dispersion re- 
lations of PbS based on an ab initio LDA electronic band 
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structure [TU] • Using these dispersion relations, which are 
rather close to the few points measured with inelastic 
neutron scattering (INS), we have calculated the density 
of two phonon states for both sums and differences of 
two phonons, which is relevant to the interpretation of 
the two-phonon Raman scattering. A thorough discus- 
sion of the observed spectra based on the calculated two 
phonon density of states is given. 



II. EXPERIMENTAL 

The samples were similar (in some cases identical) to 
those used in Ref. [8] The synthetic samples were grown 
as discussed in Refs. [HI HU The natural PbS samples 
whose spectra we report were mineral galena (cleaved 
surfaces) from Creede, Colorado. They were p-type with 
hole concentrations Nh ~ 2x 10 17 cm~ 3 . The isotopically 
modified samples were also p-type, with ~ 10 18 cm~ 3 . 
They were mounted in a closed-cycle He-cryostat (CTI- 
Cryogenics, Mansfield, MA USA) with temperature con- 
trol in the range 10 — 300 K and optical quartz windows. 
Raman measurements were performed using the 647 nm 
line of a Kr + -ion laser which was focused to a ~ 50 /im 
diameter spot with 50 mW incident power measured at 
the entrance window of the cryostat. The signal was 
collected using a photographic zoom lens (Canon, xlO 
magnification) of high numerical aperture (0.45) onto the 
entrance slit of a high-dispersion double-additive U1000 
Jobin-Yvon spectrometer coupled to a liquid N2-cooled 
CCD detector (512 pixels) without any polarization anal- 
ysis. Integration times in the range 5 — 15min (repeated 
3 times, for the software to identify and discard spuri- 
ous cosmic rays in the spectra) were used for different 
measurements in each window of the spectral range. The 
high dispersion of the double U1000 spectrometer in the 
red region (~ 647 nm) implies that several CCD- windows 
(eight) are needed to cover the full spectral range from 
150 to 550 cm -1 . Data acquisition of the intrinsically 
weak second-order spectra is therefore slow and ulti- 
mately limited by the long-term stability of the optical 
system plus cryostat. The chosen integration times were 
a compromise between signal-to-noise ratio and a total 
integration time that would allow a full temperature scan 
within a period of several hours. 

Figure [TJ shows our raw second-order Raman data for 
natural PbS at different temperatures. These spectra 
are dominated by two main features in the ~ 210 and ~ 
430 cm - regions. There is evidence in the experimental 
spectra for substructure in these peaks, to wit: (i) The 
singularity at ~ 430 cm -1 is clearly a doublet and can be 
represented as such throughout the entire temperature 
range by means of two gaussians. {ii) The singularity 
at ~ 210 cm" 1 only reveals itself as a broad doublet at 
temperatures below ~ 30 K, as we shall show in what 
follows. Above ~ 30 K it is very difficult to isolate its two 
individual components, most likely because of increasing 
anharmonic broadening [9|. 
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FIG. 1: Raw temperature dependence of the second-order 
Raman scattering of natural PbS at different temperatures in 
the range 12 — 120 K. The data are for 647 nm laser excitation. 

As already mentioned, the 647 nm Kr + -laser line was 
one of the few lines with which data could be obtained in 
both natural and isotopically substituted PbS. Going far- 
ther into the visible (in the green, for example, at 514 nm, 
see Fig. 2 of Ref. ITTjl resulted in undesirable (extrinsic) 
backgrounds that almost completely obliterated the (al- 
ready very weak) second-order spectra (see also Fig. 8 
of Ref. [T2"]l . On the other hand, going into the near-IR 
region (e.g. the 676 nm line of a Kr + -ion laser) results 
in very weak intrinsic second-order spectra. Laser exci- 
tation around 647 nm seemed to be a good compromise 
between low background and reasonable signal-to-noise 
level (possibly aided by a weak resonance effect in the 
red |llj). Another important experimental detail is the 
power level of the laser. The surface of PbS tends to 
be rather sensitive to photodecomposition (mainly, but 
not exclusively, in the presence of oxygen). The power 
densities used in these measurements are well below any 
measurable photo-degradation of the surface throughout 
the total integration time. 

The analysis presented here is based on the observa- 
tion of two main structures in the Raman spectra: (i) 
one structure at high wavenumbers, around ~ 430 cm , 
which is a doublet and has already been related to the 
generation of two optical (actually LO) phonons [H ITT] 
and, (ii) a second structure at ~ 210 cm -1 , which is pos- 
sibly also a doublet but cannot be well resolved in most 
of the measured spectra and must be related to either 
two acoustic phonons (2 LA) or two optic ones (2TO) as 
will be discussed later. Our experimental approach is to 
characterize the behavior of these spectral features as a 
function of temperature and to ascertain how they dif- 
fer between natural and isotopically substituted samples. 
Beside these main features, some weak and broad struc- 
tures have been surmised at ~ 160, 340, and 500 cm -1 . 
These structures probably correspond to those in Fig. 



3 



4 of Ref. 12 at 174, and 320 cm" 1 . In this figure, a 
rather sharp peak appears at 135 cm -1 which is below 
the frequency range of our spectra and we cannot def- 
initely assign to PbS. The peaks observed in Fig. 4 of 
Ref. [TJ]at 600 and 961 cm -1 do not seem to correspond 
to the second-order spectrum of PbS. They are probably 
due to the v\ and vibrations of SOj 2 tetrahedra in 
some form of lead sulfate [TBI [Tl] . 

Resolving the different peaks which contribute to the 
main two singularities under discussion is the first step to 
understand their temperature dependent behavior and, 
accordingly, their dependence on isotopic mass. The ba- 
sic limitation here is the lack of a reliable anaytical model 
for the singularities in the phonon densities of states (see 
next section) responsible for the second-order scattering 
observed. Gaussian lineshapes represent the individual 
singularities rather well, and will be used throughout our 
analysis. The fact that Gaussian lineshapes work well 
may also indicate that the shape of the observed struc- 
tures is heavily influenced by anharmonicity, rather than 
by the specific type of DOS critical points involved. An 
additional problem is the reliability of the fits at high 
temperatures, where even the high energy doublet around 
~ 430 cm -1 merges into a visually unresolved single peak, 
and the signal-to-noise ratio worsens considerably. Fits 
can still be obtained by restricting the range in which 
the different parameters are allowed to vary. Resolution 
of the two components, however, is not a problem at 
temperatures below ~ 60 K: at these temperatures the 
~ 430 cm -1 peaks are resolved even to the naked eye. 

We tried two ways of isolating the different contribu- 
tions to the singularities as depicted in Fig. [2] One is 
a direct fit with two Gaussians, an example of which is 
shown in Fig. [2ja) for the high-wavenumbers feature. 
An alternative way that works satisfactorily through- 
out the entire temperature range is also to estimate 
the energy of the peaks by numerically calculating the 
second derivatives of the observed spectra (a common 
practice when dealing with electronic dielectric functions 
and the related energy gaps) ; such derivative spectra are 
shown in Fig. [2Tb). In this latter case the analysis be- 
comes independent of the specific lineshape model used 
to fit the data: the critical energies can be assigned to 
the wavenumber of the deepest spectral features in the 
derivative spectra (see arrows in Fig. |2jb)). 

The best comparison between natural PbS and iso- 
topically substituted Pb 34 S is obtained from the direct 
observation of the spectra at the lowest achievable tem- 
perature in our experimental setup (i.e. 12 K), as shown 
in Fig. [3] The maximum of the low- wave numbers feature 
shifts from ~ 205 cm" 1 in natural PbS to ~ 201cm" 1 in 
Pb 34 S at 12 K (~ 2% shift). However, the correspond- 
ing lineshapes are not exactly the same and therefore 
the shift so obtained is not expected to be very pre- 
cise. The maximum of the high-wavenumbers feature, 
on the other hand, moves by ~ 15 cm" 1 , from 432 down 
to 417 cm -1 (3.5%) between the two samples. Despite dif- 
ferences in the lineshapes between one case and the other, 




Raman shift [cm" 1 ] 




150 200 450 500 

Raman shift [cm 1 ] 



FIG. 2: Analysis of the second-order Raman scattering fea- 
tures of natural PbS. Two approaches have been followed; 
both giving similar results. In (a) we show the direct fit of 
the high energy peak with two Gaussians at T = 12 K. The 
fit has a few imperfections in the lineshape but is an accept- 
able representation of the singularity in the density of states. 
The numerical second derivative approach to the estimation 
of critical point energies is shown in (b) for both singularities 
(at high and low wavenumbers) in Fig. [T] (note the break in 
the horizontal axis) . The low wavenumbers "acoustic" feature 
splits into a clear doublet at low temperatures. The energies 
of the two peaks in the high energy feature (obtained from the 
local minimum in d 2 I /du) 2 ) can be easily followed throughout 
the entire temperature range. See the text for further details. 
The arrows show the wavenumbers of the bands contribut- 
ing to the feature at ~ 430 cm" 1 at both T = 12 and 120 K, 
respectively. 



the feature at high wavenumbers seems to be more sus- 
ceptible to a mass change in S than that at low wavenum- 
bers. 

The temperature dependence of the low-wavenumbers 
feature is difficult to follow for a variety of reasons which 
include the fact that a doublet cannot be discerned above 
~ 40 K. In addition, this feature does not seem to be very 
sensitive to the sulfur mass substitution. For a discus- 
sion of the dependence of the observed shift with temper- 
ature we therefore concentrate on the high-wavenumber 
feature at ~ 430 cm" 1 . This feature shifts up in fre- 
quency with increasing temperatures (an unusual sign of 
the shift |9J) Note that the 210 cm -1 peak, and that at 
150 cm" 1 , both show the usual down-shift with increas- 
ing temperatures [5]. 
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FIG. 3: Comparison of the raw Raman second-order spectra 
between natural and isotopically substituted Pb 34 S at 12 K. 
The shift to lower energies of the singularity at ~ 430 cm -1 
can be easily seen with the naked eye. In terms of percentages 
of change in the energies, the singularity at low wavenumbers 
(~ 210cm -1 ) is less susceptible to a change into a heavier S 
mass in the isotopically susbtituted sample. 
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Data in the isotopically enriched sample have been 
taken over a smaller energy range comprising only 3 CCD 
windows with longer integration times (15min x 3), in 
order to obtain the frequency of the peaks with better 
accuracy. The measurement thus requires 45 min of inte- 
gration time for each temperature. The temperature de- 
pendence of the high-wavenumber feature in Pb 34 S can 
be seen in Fig. [4j 

Of the two components of this doublet, the one with 
the highest intensity is the easiest and most reliable to 
follow as a function of temperature. The second peak, 
at higher frequencies, becomes very weak above ~ 40 K 
and its energy cannot be very accurately determined. We 
thus concentrate on the temperature dependence of the 
former which is shown in Fig. [5] for both the natural and 
the isotopically modified samples. 

The data shown in Fig. 5 for both samples exhibit the 
characteristic temperature dependence expected from an- 
harmonic interactions (the phonon under consideration 
decaying into two virtual ones). We have used for the 
fits described below both the data displayed here and 
those reported in Ref. |5] corresponding to higher temper- 
atures (which overlap only in the ~ 77 — 120 K temper- 
ature range with our datal [Po] . The Raman shift ui as a 
function of temperature (T) has been fitted with: 
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where To is a characteristic (average) energy (in units of 
temperature) of each of the two virtual phonons respon- 
sible for the main (or average) anharmonic interaction 



FIG. 4: Temperature dependence of the high- wavenumbers 
feature in the second-order scattering of Pb 34 S. Both types 
of analysis of the spectra, i.e. fits with two Gaussians or 
numerical second-order derivatives cast similar results within 
~ 1 cm -1 for the frequencies of the singularities. The temper- 
ature dependence of the low-wavenumbers peak in this dou- 
blet is displayed in Fig. [5] and compared to the counterpart 
in natural PbS. 
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FIG. 5: Temperature dependence of the main peak in the 
high-wavenumbers feature of the second-order scattering of 
both natural and isotopically enriched Pb 34 S. The solid lines 
are fits with Eq. [I] See the text for further details. 



with each of the two "bare" frequencies adding up to the 
bare luq (i.e. the wavenumber that would be observed if 
anharmonic effects were not present) , Alu is the frequency 
renormalization at T — 0, i.e. the zero-point anharmonic 
renormalization (the real part of the corresponding self- 
energy). For T -> 0, u> -> (cj + Aw), whilst if T > T , 
lj oc T which is the classical limit. If 4Tq < luq (loq and Tq 
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Parameter 


PbS (natural) 


Pb 34 S 


uo (cm -1 ) 


428.0 


415.5 


Auj (cm -1 ) 


5.01 


4.88 


T (in cm" 1 ) 


65 


65 



TABLE I: Parameters of Eq. [T] as fitted to the data in Fig. 
[5] Typical error bars for these parameters are estimated to 
be ~ 0.5 cm -1 . See the text for further details. 

in the same units) Aw is positive (because of the sign of 
the energy denominator in the expression for the real part 
of the self-energy [9 ). Figure [5] shows the best fits to the 
data with Eq. [ljfor both natural PbS and Pb 34 S. The fit 
parameters are summarized in Table [T] These parameters 
have, actually, large error bars because the wavenumbcrs 
of the peaks change by ~ 10 cm -1 over the full temper- 
ature range and are hard to pin down to better than 
~ 1cm -1 , as can be easily surmised by visual inspec- 
tion of the raw data. Based in the information discussed 
below (vide infra) regarding the dispersion relations and 
the phonon density of states, the following picture of the 
anharmonic interaction responsible for the temperature 
shift of the ~ 430 cm -1 mode emerges. The shift is af- 
fected by virtual decay into two modes (65 cm -1 ) with 
basically Pb and very little S character. It is therefore 
reasonable to fix To to be the same for both samples, 
as was done when performing the fit. This decreases 
the number of parameters of the combined fits by one 
and produces more physically meaningful values for the 
other parameters. Our fits are compatible with a change 
of ~ 3% and ~ 2.4% in w an d Aw, respectively, with the 
sign of these shifts corresponding to ojq and Aw decreas- 
ing with increasing mass, as expected (see next section). 



III. THEORY 

The interpretation of the phenomenology just de- 
scribed is helped by calculations of the phonon disper- 
sion relations of PbS and the relevant densities of one 
and two-phonon states (the latter with k-vectors equal 
in magnitude but opposite in sign). Semiempirical cal- 
culations are available in the literature [HI [17]; they were 
performed by fitting a few parameters to a number of 
experimental frequencies obtained vs. k by inelastic neu- 
tron scattering (INS), as reported in Ref. [TTJ Here we 
use parameter-free dispersion relations derived ab initio 
using the ABINIT code[18 . Details of the calculations 
are given in Ref. 1191 which also contains the density 
of one-phonon states and its projection on the two con- 
stituent atoms of PbS. The corresponding densities of 
two-phonon states with ki = ±k 2 (sum as well as dif- 
ference of the two phonons) , were not presented in Ref. 
IT9l and are given here. The "sum" processes are re- 
sponsible for the anharmonic down-conversion into two 
phonons of lower energy [20 ; while the "difference" an- 
harmonic decay describes the process by which a phonon 
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FIG. 6: Ab initio phonon band structure calculation of PbS 
along the principal directions in the Brillouin zone (in recip- 
rocal lattice units: r.l.u). The points represent experimental 
determinations from neutron scattering jl7|. Note that some 
of the bands have a much flatter dispersion than that pre- 
dicted in the model. We shall use experimental points when- 
ever possible and only assign energies to specific bands when 
the agreement between theory and experiment is good. 



is destroyed simultaneously with a second low frequency 
phonon (which must be already in the crystal) and a 
phonon of higher energy (with respect to the initial one) 
is created. The difference processes tend to have less 
importance at low temperatures and vanish for T — > 0. 
Sum and difference processes are not only responsible for 
anharmonic phonon decay but should also determine the 
second-order Raman spectra discussed here. 

Figure [6] shows the ab initio calculated dispersion rela- 
tions of PbS along three principal directions of the Bril- 
louin zone ([100], [111] and [110]). The experimental 
points obtained with INS are represented by diamonds. 
The discrepancies which exist between experimental and 
calculated points are not surprising when one considers 
that no adjustable parameters have been used in the 
latter. The largest discrepancies are found for the TO 
bands. They may be due to anharmonic shifts which af- 
fect the experimental data but are not included in the 
purely harmonic calculations. They are in part included 
in the semiempirical calculation because of the use of ad- 
justable parameters. One can nevertheless surmise that 
the ab initio calculations should provide an overall bet- 
ter agreement to the harmonic dispersion relations than 
the semiempirical ones. Figures [7] and [8] display the sum 
and difference DOS obtained through Brillouin zone in- 
tegration from the dispersion relations of Fig. [6] The 
results in Figs. [6|8]set the framework for the following 
discussion. 
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IV. COMPARISON WITH THE EXPERIMENT 
AND DISCUSSION 

The dispersion relations in Fig. [6] show large regions 
where the TA bands are rather flat. They give rise to 
a high and rather narrow peak at ^60 cm -1 (see Fig. 2 
in Ref. IT9) |. This peak is likely to correspond to the 
frequency To=65cm -1 given in the Table. In Fig. 2 of 
Ref. [19j the partial contribution of the lead and sul- 
fur vibrations (i.e., the eigenvector components) to the 
one-phonon DOS is also given. The large mass difference 
between these two atoms produces a nearly complete fre- 
quency separation of the corresponding vibrations: up 
to ~ 120 cm -1 the Pb atomic displacements dominate 
whereas the vibrations above 120 cm -1 are basically dom- 
inated by the motion of S atoms. Hence the two phonons 
contributing to the ~ 430 cm -1 band should have a 
squared vibrational amplitude proportional to (Ms) -1 / 2 
(where Ms is the atomic mass of the corresponding sul- 
fur isotope). The change of (Ms) -1 / 2 between the two 
isotopes considered here (natural S is nearly pure 32 S) 
is 3%, a value which agrees with the fitted variation in 
Auj (Aw=2.6%). 

The densities of states for the sum and difference of 
two phonons, as obtained from the ab initio dispersion 
relations, are displayed in Figs. [7] and [8j Multiplied by 
the appropriate Bose-Einstein factors and coupling con- 
stants they determine both anharmonic self energies and 
two phonon Raman scattering intensities. The sharp two- 
phonon peak seen in Fig. [7]at ~120cm -1 corresponds to 
the two phonons (each at T =60 cm -1 ) involved in the 
virtual anharmonic decay of the 430 cm -1 Raman exci- 
tations. Several factors need to be considered when com- 
paring second-order Raman spectra with phonon sum or 
difference DOS. In general, the second-order spectra will 
resemble (in part) the sum DOS, i.e. only contributions 
of two phonons , with +k and — k wave vectors, respec- 
tively, and either the same (overtones) of different (com- 
binations) frequencies. This is, in fact, how the calcu- 
lated sum DOS was generated: it contains both overtones 
and combinations. Although there are many more combi- 
nations (5x6=30) than overtones (6) for each k, the over- 
tones seen are often enhanced by the coupling constants 
(i.e. matrix elements) . A detailed comparison of the 
Raman spectra with the two-phonon DOS is not justi- 
fied without the introduction of possibly resonant optical 
matrix elements and anharmonic self-energies. However, 
peaks and other structures observed in the two-phonon 
DOS may find their counterparts in features observed in 
the Raman spectra. The sets of phonons responsible for 
the calculated features can be identified by comparison 
with the dispersion relations. Anharmonic shifts can then 
be introduced by taking into account the difference be- 
tween calculated and INS-measured phonon frequencies 
shown in Fig. [6] 

Besides the main ~430cm -1 band, there is also in 
the measured spectra the already mentioned ~210cm -1 



band, (note that its frequency decreases with increasing 
T's[8], i.e. the usual behavior for the shifts of phonons 
with temperature) and a few other weak bands already 
reported in Refs. [El QUITS]- 
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FIG. 7: Phonon sum density of states (DOS) calculated from 
the dispersion relations in Fig. [6] High wavenumber features 
related to the LO optical branches have predominantly "S" 
character, while low wavenumber "acoustic-like" features have 
mainly "Pb" character. See the text for further details. 




FIG. 8: Two phonon difference density of states (DOS) calcu- 
lated from the bands in Fig. [6] The difference spectra should 
occur at the limit of what we can measure at low frequen- 
cies in our experimental setup and should also have a distinct 
temperature dependence disappearing for T —> (which is 
not observed experimentally on the low energy limit of our 
observation window) . The experimental evidence available at 
this stage suggests that the observed spectra are dominated 
by two-phonon sum processes. 

We make the following tentative assignments of both 
the main and secondary features: 

• The experimental curves have a small (but dis- 
cernible) feature at about ~ 160 cm -1 . This is at 
the low energy limit of our experimental range and 
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therefore can only be assigned tentatively. It could 
correspond to TA + TO (combination band) along 
E (~ (55+ 100) cm" 1 = 155 cm- 1 ). The frequency 
of this spectral feature also decreases with increas- 
ing temperatures [5]. 

• The experimental peak in the region 205—210 cm -1 
seems to have a shoulder at ~ 230 cm -1 seen clearly 
at low temperatures. This feature was previously 
attributed to LO scattering (which is formally for- 
bidden, but can be made allowed and observable 
by Frohlich interaction when the laser frequency 
is near an interband resonance [TT]). The data in 
Refs. 0E], which were not taken at low temper- 
atures, did not allow a clear identification of the 
shoulder at about ~ 220 cm -1 which could be due 
to two-phonon scattering. In fact, both features at 
~ 205 — 210 cm -1 could correspond to two-phonon 
scattering, and it is hard to distinguish this assign- 
ment from that to a forbidden LO scattering situ- 
ation. The sharp peak at ~ 200 cm -1 in the calcu- 
lated sum DOS corresponds to 2TO phonons either 
along A or E. One could also take along E the ex- 
perimental TO+LA values and obtain ^190 cm . 
In addition, 2TO at both A and E should give a 
feature at ~ 200 cm -1 . With suitable weighting by 
matrix elements, all these assignments may play 
a role in the Raman feature observed around 210 
cm - 1 . 

• Combination scattering of LA + LO along E would 
give (from the experimental values) ~ (100 + 
230) cm" 1 = 330 cm- 1 , or ~ 340 cm- 1 along E. 
This phonon combination could be responsible for 
the observed experimental weak structure at ~ 
350 cm- 1 . 

• The broad peak in the calculated sum DOS cen- 
tered at ~ 270 cm- 1 could be twice the (theoret- 
ical) low TO along sigma and have no connection 
with the experimental phonons (there are no INS 
points around ~ 135 cm- 1 and therefore we do not 
expect a feature at 270cm- 1 in the experiment). 
The actual TO seems to be lower, possibly because 
of the anharmonic shift. 

• Finally, the asymmetric broad experimental band 
at ~ 430 cm- 1 , whose temperature dependence we 
discussed in the previous section (Fig. [5|, could 
correspond to (using experimental values (INS) for 
the phonon energies of the bands) : (i) 2LO at A = 
2 x 230 cm- 1 = 460 cm- 1 ; (it) 2LO at A = 2 x 
245cm- 1 = 490cm- 1 ; and/or (in) LO + TO at L 
and A = (190 + 240) cm" 1 = 430 cm" 1 

• For the two-phonon difference DOS in Fig. [8] there 
is less experimental evidence. Unlike the sum DOS 
of Fig. |7J the intensity of the features in the dif- 
ference spectra should vanish at low temperature. 
Any peak whose intensity vanishes as T approaches 



the lowest possible temperature (12 K in our case) 
is a strong candidate to be compared with the dif- 
ference DOS. We do not see any feature that could 
safely be assigned to the difference DOS, which ex- 
hibits strong structures at 30, 90 and 140 cm- 1 . 
All these structures lie below the lower limit of the 
spectral range that can be recorded with our CCD 
detector without being "swamped" by stray light 
from the laser. The sharp peak in Fig. 4 of Ref. 
[12] could be assigned to the peak seen at the same 
frequency in Fig. [8] (due to LA — TA phonons at 
X and A) if it is confirmed that it belongs to PbS. 
The best way to buttress this conjecture would be 
to measure its intensity as the temperature is low- 
ered. Unfortunately, Ref. jT^j only reported data 
obtained at room temperature and our spectrome- 
ter does not operate below ~ 135 cm- 1 . 

• Another aspect we have already touched upon con- 
cerns the interpretation of the unusual frequency 
up-shift of the ~ 430 cm -1 feature in Fig. |5]due to 
virtual anharmonic decay into two phonons. The 
sign of this coefficient, corresponding to a "blue" 
shift, is related to the fact that the two-phonon in- 
termediate states are at lower frequency than each 
of the phonons under consideration (2 x 60 cm -1 = 
120cm- 1 vs. ~ 215cm- 1 . The ~ 215cm -1 modes 
are all S-like, whereas the 60 cm -1 phonons are 
lead-like. It has been mentioned in the previous 
section that this leads to a zero-point anharmonic 

— 1/2 

renormalization Aw proportional to M s . The 
corresponding renormalizations for the two isotopes 
measured are given in the Table. They differ by 
2.6%, a value which agrees with that determined 
from the corresponding isotopic masses of sulfur. 
For the 205 cm- 1 peak, regardless of whether it is 
due only to forbidden one-LO scattering or to two 
phonon scattering, each with a frequency around 
100 cm- 1 , most of the spectral weight of the two 
phonons with which it interacts lies above the fre- 
quency of the phonons being considered (see Fig. 
0. There is some experimental evidence for this 
down-shift in our data, (see also Ref. [TT]) even 
though the resolution is poor and the exact nature 
of the ~ 210 cm- 1 peak in terms of one-phonon for- 
bidden scattering or two phonon scattering cannot 
be unambiguously decided at this stage. Our anal- 
ysis, however, paves the way for the interpretation 
of future experimental observations. 

• The two components of the ~ 430 cm- 1 band must 
be sulfur like, according to Figs. [6] and [7J Hence 
they should shift (according to the M -1 ^law) by 
3% from natural PbS to Pb 34 S. The measured shift, 
as listed in the Table, is also 3%. We note that we 
have estimated a somewhat smaller shift (2%) for 
the measured spectra of the ~ 210 cm- 1 band. Al- 
though the difference between this shift and that of 
the square roots of the sulfur masses (3%) lies prob- 
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ably within the experimental error, it would be easy 
to justify it in terms of the theoretical expectation: 
For the one-LO scattering the shift should be deter- 
mined by the reduced mass of Pb and S and amount 
to 2.6%. In the case of scattering by two phonons, 
with frequencies around ~ 100 cm -1 , which fall in 
the region which separates Pb-like from S-like vi- 
brations, a shift smaller than that predicted from 
the sulfur mass can be expected but it is difficult 
to estimate quantitatively. 



corresponding densities of two-phonon states. The main 
spectral feature, abroad doublet centered at ~ 430 cm -1 , 
has been assigned to scattering by two LO phonons. 
Its "anomalous" temperature dependence has been at- 
tributed to virtual anharmonic decay into two phonons 
of frequency ~ 65 cm -1 each. The dependence of the fre- 
quencies of the two main observed features on the isotopic 
mass of sulphur has been shown to agree with theoret- 
ical predictions, a fact which we believe buttresses our 
assignments and interpretation. 



V. CONCLUSIONS 

We have measured the Raman spectra of two PbS sam- 
ples, one made out of natural galena and the other a syn- 
thetic crystal with natural sulfur replaced by S, in the 
150-500 cm -1 spectral region at temperatures between 10 
and 120 K. These data extend earlier data on the same 
materials covering the range 100-300 K. For the quanti- 
tative interpretation of the spectra we have performed ab 
initio calculations of phonon dispersion relations and the 
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